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Abstract 

A complete list of operators contributing at the lowest order to Quartic Neutral Gauge 
Boson Couplings involving photons and Z-bosons, is presented. We show that, for the 
couplings we consider, the lowest order contribution is from dimension 8 operators in 
the case when a light Higgs is present and from dimension 6 operators in the higgsless 
case where electrowcak symmetry is non-lincarly realized. We also show that these 
operators are generated by exchange of the Kaluza-Klein partners of the graviton in 
extra-dimensional models. We then explore the possibility of probing these couplings in 
the diffractivc photon fusion processes pp(77 — >■ 77)pp and ^^(77 — )■ ZZ)pp at the 14 
TeV LHC. We find that the 7777 coupling can be probed most sensitively and values 
as small as 1/(1.8 TeV)^ can be measured. For the 77ZZ coupling, values as small as 
1/(850 GeV)^ and 1/(1.9 TeV)^ can be probed in the light Higgs and higgsless cases 
respectively, which is an improvement by orders of magnitude over existing limits. 
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1 Introduction 



The Standard Model (SM) has been tested very accurately by experiments. There are, how- 
ever, many theoretical reasons to believe that there is physics beyond the SM. Some of these 
motivations, like the hierarchy problem and the existence of dark matter, point to the exis- 
tence of new physics at the TeV scale. LEP-2 precision data and flavor constraints seem to 
favor a scenario with a mass gap between a light Higgs {rriH ^ 200 GeV) and new physics at 
the scale of a few TeV. A model independent way of parametrizing the effects of new physics 
in such a scenario is to use the effective field theory approach. All possible operators allowed 
by the symmetries of the theory are included, suppressed by appropriate powers of the cut-off 
A. If A is the order of a few TeV, these operators can be directly measured at the Large 
Hadron CoUider (LHC). 

These operators are expected to give rise to anomalous triple (T||2] and quartic gauge 
couplings [3]. In this work we discuss a special class of these couplings: the Quartic Neutral 
Gauge Boson Couplings (QNGC), that is, quartic vertices involving only the neutral gauge 
bosons, 7 and Z. QNGCs are special because as we will show they do not exist in the SM and 
receive their lowest order contributions from dimension 8 operators. Thus the measurement 
of these couplings would indicate directly the presence of dimension 8 operator^ in a scenario 
where a light Higgs is present. This is not true, for example, in the case of 'j'jW~^W~ and 
ZZW^W~ couplings which get contributions from the SM lagrangian and its dimension 6 
extension [l]. Thus QNGCs can be very useful in probing new physics scenarios with a light 
Higgs that exclusively generate dimension 8 operators. One such example that generates only 
dimension 8 operators at tree-level is the exchange of the spin 2 Kaluza-Klein excitations of 
the graviton in models with large extra dimensions. We will see how integrating out these 
massive modes generates QNGCs and how probing these couplings would allow us to probe 
the fundamental Planck scale in these extra dimensional theories. 

We also consider the higgsless case where electroweak symmetry is non-linearly realized. In 
this case, with the exception of the ZZZZ coupling, QNGCs do not appear at the dimension 
4 level and the lowest order contribution comes from dimension 6 operators. This is unlike 
quartic gauge couplings having W^^^ bosons which always appear first at the dimension 4 
level. Thus in this case also, unlike processes involving quartic gauge couplings with W^^~ 
bosons, processes involving QNGCs can directly probe higher order operators (in this case 

^In the case of the couphng there is a non-local contribution from the 77 — > /i* — > ZZ process which 

is of an order lower than dimension 8 contributions. However, as we discuss later, this contribution can be 
subtracted if the h — 77 partial width is known. 



2 



Figure 1: The diffractive photon fusion processes ppi'yy ll)pp and ppi'yy — )■ ZZ)pp. The 
outgoing protons can be detected by very forward detectors to be installed by ATLAS and 
CMS. In the figure above Oi represents operators contributing to Quartic Neutral Gauge 
Couplings (QNGCs). 

dimension 6 operators). 

In this work we will explore the possibility of measuring the 7777 and 77ZZ couplings 
in the diffractive photon fusion processes, pp(77 — )■ 77)pp and pp(77 — )■ ZZ)pp (see Fig. [T]) 
respectively. There are plans to install very forward detectors by the ATLAS and CMS 
collaborations ji] which can detect protons that scatter diffractively at small angles and thus 
can identify such processes. To the best of our knowledge, this is the first study on the LHC 
sensitivity of the measurement of the 7777 coupling. There have been previous studies for the 
77ZZ coupling, but these studies focussed only on the higgsless case. In Refs. |5|[6] probing 
the 77ZZ couphng by inelastic processes like pp — )■ 777, pp — )■ 77^' and pp — )■ jjiZZ — 
77) — )■ jj77 has been studied, whereas, in Refs. [7|[8] measurement of this coupling in the 
diffractive process ppicil — > ZZ)pp, that we will study in this work too, has been explored. 
In this work, however, we consider both the light Higgs and the higgsless cases. As we will 
see, in the higgsless case considered in the previous studies only a subset of all the operators 
relevant to the light Higgs case are important. 

Let us now see what are the advantages of diffractive photon fusion processes in measuring 
these couphngs. Even if a process can be traced back to a definite set of operators as is the case 
here, it is rarely the case that a particular collider signature can be traced back to a unique 
process. For this reason many different, complementary measurements are usually required to 
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uncover the underlying new physics processes. For example consider the inelastic counterpart 
of the signature we are considering for the coupling, the pp — ?■ jj (77 — )■ ZZ) — t- jjZZ 

process or the similar vector boson fusion (VBF) process pp — ?■ jj{ZZ — t- 77) — )■ jj'j'j. 
Although these signatures would have a much larger cross-section than the diffractive signature 
we are considering, if an excess is observed in the jjZZ or jj77 final states it would be 
hard to reconstruct the exact process responsible for it because of the many different new 
physics processes in addition to QNGCs that can have this signatur^ The triple gauge boson 
production processes pp — )■ 77Z (studied previously in Ref. [5||6]) and pp — )■ 777 are somewhat 
better in this respect but, again, because the intermediate state in pp — )■ 77^(777) cannot be 
known, it would not be possible to conclude with certainty that QNGCs are responsible, if an 
excess is seen. The diffractive signals we study in this work are interesting because exclusive 
final states where two protons have been detected in the forward detectors can arise only 
from diffractive photon fusion or exclusive pomeron fusion. In the latter case the underlying 
subprocess would be gg — '-ff/ZZ. Thus the inverse problem of pinning down the new physics 
responsible for an excess in the ppy'y and ppZZ final states is relatively less ambiguous as 
there are only two new physics possibilities namely the enhancement of the 77 — )■ 77/ZZ 
processes and/or the enhancement of the gg — )■ jj/ZZ processes. As we will discuss later, 
exclusive pomeron fusion processes are, however, expected to have a much smaller cross-section 
as compared to photon fusion processes. 



2 Operators that give rise to Quartic Neutral Gauge 
Boson Couplings 

2.1 Light Higgs case 

We want to write down the lowest order contribution form higher dimensional operators to 
QNGCs, that is, quartic vertices involving only the neutral gauge bosons 7 and Z. We will 
consider only CP conserving operators here hence we will not use the dual field strength 
tensors like 5^,^. In gauge invariant operators the electroweak vector boson fields can appear 
either from the field strengths -B^jy and W^^ or covariant derivatives acting on the Higgs field, 

■^Thc pp — !■ — > ZZ) jjZZ process is experimentally challenging for a separate reason too which 
is that the two jets would not have the special properties of VBF jets. VBF jets have a large rapidity gap 
between them and have high pt (see for example pg. 1271-1305 of [29J). For reasons mentioned in Section [5] 
the Pt of the jets in the process pp — — >■ ZZ) jjZZ is approximately equal to the photon virtuality 
and thus expected to be very small. So while the photon fusion jets would have a large rapidity gap too, they 
will have very low pT 
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that is factors like, 

D,<^ = d,<^ - ig^Wl^ - t^B,^, (1) 

provided the Higgs field, $, gets a vacuum expectation value (VEV). In the above equation 
T/ are the Pauli matrices. Note that the ejjxWjWK part of Wj^^ cannot contribute to quartic 
neutral gauge boson couplings as we must have / = J = 3 to get photons or Z fields unlike 
for instance in the case of the 'j'jW~^W~ coupling, where two gauge boson fields can come 
from the same field strength tensor VF^^. This has the important implication that the lowest 
order contribution to quartic neutral gauge boson couplings comes from dimension 8 opera- 
tor^ because we need either a field strength or -D^$ factor, both dimension 2 operators, for 
each of the four gauge bosons. There must be an even number of -D^$ factors along with 
the field strength tensors in these operators because otherwise there are an odd number of 
Lorentz indices in total and it is impossible to contract all of them. Thus we see that the 
operators we are interested in can have either four covariant derivatives, two field strengths 
and two covariant derivatives or four field strengths. We thus obtain the following lagrangian 
of dimension 8 operators. 




Note that, in the above hst, the operators BpyeuxW^P^'W-l^W^P'' , Bp^eijK{^^a^^)W-l^W'^P'' 
and B p^eij kW^ PPW •^^"W are absent because they are all equal to zero. In the first two cases 
eijK is antisymmetric in J and K whereas the rest of the operator is symmetric in J and K 
and in the third case Bp^ is antisymmetric in /i and v whereas the rest of the operator is 
symmetric in and v. Operators with two a^s do not appear above as these can be reduced 
to operators in our list using the identity, a^cy^ = 5^^ + ie^^^ax- Also notice that there 
are no operators like ^WPD^^BppB^P . This is because such operators can be expressed as 
linear combinations of total derivatives, operators already in the list and operators that do 

''This fact is also true for 3577 and ggZZ couplings (g being the gluon). 



5 



not contribute to QNGCs, as follows, 



-^^D^^d^iB^pB'P) (3) 

where we have used dp_{^^Dj,^) = D^^W^^ + ^^{D^D^^). Finally, operators with two 
Levi Civita tensors, like ep^ypaeap-ysBP'^B^^ B^^ B"^ , which corresponds to taking two dual field 
strength tensors, are not in the list. This is because using the identity, 

we can express such operators in terms of operators contracted using metric tensors which are 
already in our list. 

If we rewrite these operators in terms of the fields A and Z defined by. 



and the field strengths. 



Ws [ [ Z 



Fpu — dpAu — dyAp 



(5) 



Zpy = dpZy - dyZp. (6) 

we will get ZZZZ, ^ZZZ^ ^^ZZ ^ 777Z and 7777 couplings. In this work we will explore 
the possibility of measuring these vertices by diffractive processes involving photon exchange, 
that is processes like pp(77 — )■ X)pp. Such processes can measure the •yyZZ, 777Z and 7777 
vertices but we will focus only on the 77ZZ and 7777 vertices here (we give the complete list 
of the 777^', 'yZZZ and ZZZZ couplings in Appendix A). Expressing the operators above 
in terms of A and Z fields and the respective field strengths using eq. |5| [6]we get, 

77 

ClVrlr^ = ^F„,,.FP''Fp^FP'' + ^Fp,F''PFp^F'''' 



'7777 _ 
■QNGC — 


77 

^F 1 






■QNGC ~ 


A4 2 

zz 

+ ^ F 




' A4 



zz /i#2 ZZ 

CV'nn = "^'-^F.^FP'^ZpZP + j^^Fp^F^^PZpZ^ + j^Fp^FP^Zp^ZP^ 



zz zz 
p.ZP-Fp^ZP'^ + -j^Fp.FPPZp^Z'^^ + -j^Fp^ZPPFp^Z'^'' (7) 
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where, 





= c^cs + s^cq + clsKdo + Cii) 




oT 


= c^cis + stci4 + C^S^(ci5 + Cie) 






_ 2 _|_ 2 




ar 


2 1 2 
— C^Ce + S^Cf 






= 2C^S^(C8 + Cg + Cio) + (st + ct)cii 


^ 2c^S^Ci2 


ar 


= {cl- slfci2 + 4:clsl{c8 + Cg + Cio) 


- 4c^s^cii 


ar 


= 4:clsl{ci3 + ci4 + cis) + (st + ct - 


2c^4)ci6 - 


ar 


= (ct + ■st)ci7 + 2c^S^(ci3 + Cl4 + Ci5 


) - 24s^ci6 



We have thus hsted all operators that contribute to the 7777 and 77ZZ' vertices. Note that 
Ci and C2 do not appear in the RHS in eq. |8] because the corresponding operators contribute 
only to the ZZZZ coupling. As we want to measure these couplings by the 77 — )■ 77 and 
77 — )■ ZZ processes respectively, let us also list operators that might enhance non-local 
background contributions through the processes like 77 — ?■ X* — )■ 77/ZZ at the same order, 
where X is some SM field. We find that the only dimension 6 operators giving such non-local 
contributions at the same order are those that contribute via the 77 — > /i* — ZZ processes due 
to the anomalous h'yy couplings they introduce. These operators (already listed in Ref. [l]) 
are, 

Ce = ^m'^B.^B'^'^ + ^Ic^l^wl^^w"-" + ^i^^ai^)B,,W'^''. (9) 
which give the following dimension 5 operator once the Higgs field, $, gets a VEV, 

"^hF,.F^^ (10) 

where v = 246 GeV and, 

ao = clbi + slb2 - CwSnjbs- (11) 

Note that heavy particles that interact with the Higgs boson and photons generally induce 
this operator when integrated out. Fortunately, the coupling oq can be accurately measured 
at the LHC by measuring the /;,—)■ 77 partial width. Thus the effect of the only dimension 6 
operator that contributes to the cross section can be subtracted. 

We have not identified couplings for 77X or XZZ, that arise from dimension > 6 operators 
and contribute to this process by an X exchange, because these contributions would have the 
dependance ~ l/A*^ with n > A which would be of higher order than dimension 8 which is the 
lowest order at which the QNGCs get a contribution. 
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2.2 Higgsless case 



Note that our treatment in the previous subsection differs from that in Refs. [5j-|8] where 
only the 77 — i- ZZ process has been discussed. Among all the terms in eq. [7] the authors 
consider only the operators (9f^ and Of'^ (i.e., the operators that have the coefficients af^ 
and af ^ respectively). This can be justified if there is no light Higgs and electroweak symmetry 
breaking (EWSB) is non-linearly realized at low energies. Let us see why this is so. We follow 
the construction of Burgess et al. (Ref. lo]), use the matrix, 



E = exp(2iXi7r,/w) (12) 

and the covariant derivative, 

D;L = S^D^S - tE^[gW;Ta + g'B^Y]J:. (13) 
to define the following fields, 

eA^ = 2«Tr[Xe„D^S] 

gW^ = 2v^Tr[T±D^S]. (14) 

Here Y is the hypercharge generator, T± = Ti ±T2, where Tj are the SU(2)i generators. Xem 
and X3 are orthogonal linear combinations of Y and T3, being the unbroken generator 
of U(l)em- We have kept the unconventional normalization of Ref. |9^|, viz Tr[Ta,Tb] = ^6ab, 
TT[Ta,Y] = and Tt[Y^] = \. 

As shown in Ref. |9], the fields Z and in eq. 14 transform purely electromagnetically 



and exactly like A, Z and respectively. In the unitary gauge, S — 1 so that A ^ 
Z ^ Z and — ?■ . It is thus easy to construct gauge invariant operators we are 
interested in if EWSB is non linearly realized. In the unitary gauge these are just all possible 
operators constructed from the A, Z and fields that respect the U(l)em symmetry. We 
get therefore for the ^^ZZ coupling, 

r.HiggsleSS (fi'/2Cu.) Cll „ jpUiy ryp I (5'/2Cio) a2 T?llpry ryv /irN 

l-QNGC = -j^ ^M^^ ^P^ -j^ ^/^^^ ^P^ l^^J 

Note that in this case we get as the lowest order contributions to QNGCs two dimension 6 
operators, which are same as as Of^ and Of ^ in eq. [7| and none of the other operators in 
eq. [7] are present above. Thus, unlike in the case with the light Higgs boson, these operators 
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are indeed more important here, and this is why they are the only ones that appear in the 
analyses of Refs. [5]-[8]. The 7777 coupling does not get any contribution at this order. 

Another way to understand the above fact is by using the goldstone boson equivalence 
theorem which states that at high energies longitudinal gauge boson production processes 
should have the same amplitude as processes in which the corresponding goldstone bosons are 



produced. The operators in eq. 15 arise from operators like c Tr[{D^J:y D^'J:]F^^F'"' . This 
can be expanded to give the terms involving the goldstones like, {{c/ A'^){dpTfd^7r/v'^))Ff^^Ff^'^ , 
which tells us that the ^(77 — )■ ZlZl) amplitude will be (9(cs^/(f ^A^)), ignoring the di- 
mensionless electroweak couplings. This is larger than the amplitude due to an operator like 
((c/A^)Zp^Z'"^)F^^F^^ not included in eq. M which will give ^(77 ^ ZZ) = 0{cs^/A'^). 



The crucial difference is that, unlike the light Higgs case, the goldstones here are strongly 
coupled and suppressed by factors of 1/f and not 1/A. 

2.3 Graviton exchange in extra-dimensional theories as a source of 
QNGCs 

In extra-dimensional theories where the fundamental gravity scale can be a few TeV, the 
graviton is accompanied by Kaluza-Klein (KK) partners in the 4D effective theory. Exchange 
of the (4 + 5)-dimensional graviton, 6 being the number of extra dimensions, can be thought 
of as the excahnge of the 4D graviton and its massive KK partners. The effective operator 



induced by tree-level graviton exchange is given by 10 



47r /T T^"' 1 T^T^'X 

where T^'^ is the energy-momentum tensor. At tree level only dimension 8 operators are 
induced (at loop level only one dimension 6 operator operator is induced by virtual graviton 
exchange but this is a four fermion operator not involving the gauge bosons or the Higgs flS]). 

Almost all the operators in eq. [2] can be obtained by expanding T^'^ in eq. 16 To show 



this let us write down the energy-momentum tensor for the i?^ and gauge bosons and the 
Higgs boson. 
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Note that virtual graviton exchange will also generate operators involving the gluon field 
strength, Gj^^, like G^^^G^^"" Bp^Bf" , Gj^^G^^"' Dp^W^^ etc. Such operators would enhance 
the signal by contributing to the central exclusive pomeron fusion process (CEP), pp{CEP — )■ 
'j'y /ZZ)pp. The luminosity of photons produced by the protons is however higher than the 
luminosity of the pomerons produced that undergo exclusive fusion (by exclusive we mean 
that the pomerons do not disintegrate into fragments) by a few orders of magnitude at the 
high energies where these operators become important (see Fig. 2 in Ref. [ll]). Thus the 
ppiCEP — )■ '~^'~^/ZZ)pp contribution is expected to be negligible compared to the pp(77 — )■ 
/ZZ)pp contribution. In any case any contribution form this channel would only enhance 
the signal and thus improve the experimental potential of observing effects of virtual graviton 
exchange. 



3 Constraints 

QNGCs are very weakly constrained by existing data. There are no constraints on 7777 
couplings and the only constraints are on 77ZZ couplings. We first consider the light Higgs 



case discussed in Section 2.1 A LEP analysis |12| based on the e'^e — > Z77 process puts the 



following constraints on the operators Oi and O2 in eq. [7| 



1 aP 1 



(69 GeV)4 A4 (93 GeV)^ 
and. 



- < ^ < - (19) 

(65GeV)4 A4 (eSGeV)^' ^ ' 

While the authors of Ref. [12] did not carry out their analysis for the other operators in 
eq. [7| as these are also dimension 8 operators we expect their contribution to these processes 
to be of a similar magnitude. Thus the constraints on these couplings are also expected 
to be very weak. Somewhat stronger constraints can be derived from electroweak precision 
data. In Ref. (6] precision constraints on the operators are derived and they find the bounds, 
|aff/A^| < 1/(270 GeV)^^, which, as we shall see later, are still far too weak compared to 
the expected LHC sensitivity. 



zz 



As the operators in the higgsless case discussed in Section 2.2 are exactly the two operators 



discussed above, the only difference being that we use a different parametrization for the 



couplings, the same constraints can be translated to the couplings in eq. 15 in the higgsless 
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case, 



and, 



1 



(27GeV)2 
1 



< ^ < 



,hl 



(50 GeV)2 



1 



(24GeV)2 ^ A2 ^ (24GeV)2' 



(20) 



(21) 



whereas the precision constraints in Ref. |6] imply \a^yA^\ < 1/(420 GeV)-^. 

Now we discuss the constraints on the scale, At, for virtual graviton exchange, which 



appears in eq. 16 The strongest constraints on A^^ come from LHC data at 7 TeV. With 36 
pb^^ CMS data at 7 TeV the pp — t- jj process can be used to derive the constraint A^ > 3.8 



TeV |[14| at 95 % confidence level. The same process puts the constraint A^ > 3.6 TeV 14 
with 36 pb"^ ATLAS data at 95 % confidence level. With 1.1 fb"^ CMS data, the pp 77 



process puts the weaker constraint A^ > 3.1 TeV at 95 % confidence level [15j, but this process 
is eventually expected to probe scales up to about A^ = 6 TeV 10 . 



Finally, consider the operators in eq. [9j At tree level the couplings bi and 62 in eq. |9] 
renormalize the coefficient of the kinetic terms for the gauge bosons and W^^ which is 
equivalent to a renormalization of the couplings g' and g. Thus all tree level effects due to 
the bi and 62 can be absorbed in a redefinition of the couplings and hence these couplings are 



unconstrained. The coupling 63 is related to the S'-parameter by 16 

rA5 



h 

A2 



Op 



(22) 



4:Sij[jC^V 

Here aem is the fine structure constant v = 246 GeV and Sw, are the sine and cosine of the 
weak mixing angle. The bound on the S'-parameter for rrih = 113 GeV and with no restrictions 



on the T parameter is \ AS\ < 0.3 16 at 90% confidence level. This translates to the following 
bound on 63, 

h 



A2 



< 



1 

(6.6TeV)2' 



(23) 



This coupling is also constrained by measurements of the triple gauge couplings but these 



constraints are far weaker 17 



4 High energy behavior of amphtudes and violation of 
unitarity at tree-level 

First let us look at the 77 — > 77 process. We can find out the high energy behavior by 
dimensional analysis. The high energy behavior of the contribution from the local operators 
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in eq. [7] differs from the contribution from tlie non-local process, where the /177 vertex is 
derived from the operator Oq, as follows, 

S2 



VS 1 VS 2^^-^ 



Co : A{-f-f ^ h* ^ -f-f) ao — - X ao— aQ — 

(24) 

where s is the photon-photon center of mass energy squared. The local contribution is thus 
expected to dominate over the non-local contribution at high energies. 

For the 77 — )■ ZZ process the operators in the light Higgs case in eq. [7] can be divided into 
into two categories according to the final polarization of the Zs. At high energies the operators 
Of^, of contribute mainly to the production of longitudinally polarized Z-bosons through 
the process 77 ZlZl while the operators O^^ — O^^ contribute mainly to transverse Z 
production through the process 77 — i- Z^Zt- This can be understood by using the goldstone 
boson equivalence theorem. The operators like Of ^ — Of ^ do not arise from dimension 8 
operators involving the Higgs field (see eq. |8]) and so they do not introduce new couplings 
to the Goldstone bosons (that are eaten by the gauge bosons in the unitary gauge). New 
contributions to the process 77 — t- Z^Z^ they introduce are, therefore, suppressed. 

Let us now see the energy dependance of the 77 ZZ amplitude of the dominant Z- 
polarization modes for the different operators using dimensional analysis in the high energy 
limit, 

Ol^-Ol^ : A{^^ ^ ZtZt) ^ a~ 
Ol^Or : ^(77^^L^L)~a.^^~a.|^ 

Co : ^(77-^^ ~^ ^tZt) ^ ao — ^— 



2 



Oo : ^(77^/^ ^^^^^)~"°A^J^M|~«°A^- 

(25) 

where the the s/M| factor for the longitudinal modes comes from the longitudinal polarization 



vectors and gMz/cu, is the SM hZZ couphng. Note that according to eq. 25 the 77 — )■ /i* — )■ 
ZZ process would mainly produce longitudinal Zs. As discussed earlier the contribution of 
the operator Oq can be subtracted by measuring the /;.—)■ 77 partial decay width. For the 



operators in the higgsless case in Section 2.2, the dominant mode will be 77 — )■ Z^Zi and the 
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energy dependence would be, 

^(77 ^ Z,Z,) ~ (^7/2c.)^af ^ (26) 

z 

As all the amplitudes above grow with energy they would all violate partial wave unitarity 
for some value of s. We obtain the perturbative unitarity bound for the processes in Appendix 
B. The condition that perturbative unitarity is not violated is, 

(Re(Mf + /3^|a,r + 5,<^, (27) 



where a; (bi) is the l-th. partial wave amplitude for the 77 — )■ ^^(77) process, /3 = y 1 — s 
is the photon-photon center of mass energy, 61 is the positive contribution from other processes 
and €3 and €4 are the polarizations of the Z bosons produced. For the first term the final 
polarizations are same as the initial. The initial polarizations of the photons must be chosen 
to maximize the LHS to get the most stringent possible bound. We find the most stringent 
bounds from the / = mode. 

To ameliorate the growth of the amplitude with energy we can use form factors as follows. 



1 

1 + (JjAjy 



A^A{ ^ , ,,,,,, I . (2^ 



By Taylor expanding the modified amplitude we see that by introducing a form factor we 
effectively introduce higher order contributions, such as those expected from loop effects and 
higher dimensional operators, to cancel the tree-level growth of the amplitude. For example 
if ^ = ks'^/A'^ the for the choice m = 2 and n = 1 the modified amplitude is , 

ks^A^ 



, = (ksVA^) (1 - + {s/A)Y - [s/A)f..) . (29) 



In order that the eq. 27 is obeyed we would require that. 



{Re{hi)f < 0.1 (30) 
|ai|^ < 0.1. (31) 



63, S4 



The RHS in the two equations above do not add up to the RHS of eq. 27 because we have 
made some allowance for other contributions to 5i. To ensure that these conditions are obeyed 
we use form factors for both the 77 — )■ 77 and the 77 — ?■ ZZ amplitudes. 
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Figure 2: Here we plot (^{Re{bi)f + PYl 
mass energy for I 
Ci/A^ = 1 TeV"! We take A 
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m eq. 



vs the photon-photon center of 
with and without the form factor in eq. 28 taking 6j/A^ = 1 TeV~^ and 
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than the lower solid line where we take Ay^^^ 
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Fig. 2j shows the growth of (^(Re(6;))^ + /5 k/Pj in eq. 27 for / = with energy, for 

different choices of the form factor parameters. We consider the ligEt Higgs case in Section 2.1 
taking the couplings foj/A^ = 1 TeV"^ and Ci/A^ = 1 TeV"^. Let AJJ*^ and A^| be the values 
of s where the conditions in eq. 30 and 31 are respectively violated when no form factor is 



applied. We can see that the amplitude keeps growing for m 



n 



1 and A 



n/zz 

f 



A 



-Y-t/ZZ 
UB ' 



thus violating the perturbative unitarity bound. However, the amplitude is suppressed below 
the bound for m = 1, n = 2 and m = 2, n = 1 for the same values of AP''^^. We see that in the 



latter case the amplitude saturates the bounds in eq. |30]and[3T]at high energies. We also show 
a curve with m = 2, n = 1 but Ay/^^ = 0.9 A^l^^ which coincides with the A]^/^^ = A^/,^^ 
curve at low energies but deviates from it for s close to Afj^. Unless otherwise mentioned 
from now on we will use form factors with m = 2,n = 1 and AZ"'^^^ = AV/,^^ . While our final 



results will depend on this specific choice of form factor a different form factor would result in 
a cross-section with a different numerical value but the same order of magnitude. Thus there 
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Figure 3: The Luminosity function dL/dW in eq. 36 taking q^^x = 2 GeV^. 



will be a relatively small difference in our final sensitivity results on A as the cross-section 
goes as 0" ~ A"^ (cr ~ A"^) for the dimension 8 (dimension 6) operators in the light Higgs 
(higgsless) case. 



The Equivalent Photon Approximation and the pro- 
ton level cross-section for pp{^^ — > ZZ)pp 



Consider a general process ppi'yy X)pp. To find the amplitude for this process we need to 



know the ppj vertex. From gauge invariance the most general form for this vertex is 18 



— ^e 



(32) 



Here efip/2mp is the proton magnetic moment with /i^ = 7.78, qf, the invariant mass of one 
of the photons is, as shown in Appendix C, always space-like in such a process and thus it is 
negative. The functions Fi{q'^) and F2{q'^) in the vertex can be expressed in terms of the the 
empirically determined elastic electric and magnetic form factors for the proton, GE{q^) and 
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Ghii'f') respectively, as follows, 

Ge 



Ge + tGm 
1 + r 

Gm/ fip 



Gm — Ge 
fi:(l + r) 



:i - q'/q^,)-\ (33) 

Here GE^q^) and Guiq^) have been written in the dipole approximation with = 0.71 GeV^ 
and T = (— g^)/4mp. By a change of variables the final phase space integration for the process 
PPill ^ X)pp can be done over d{—q\)d{—ql)dujiduj2, instead of the usual variables 19 , uji 



being the energy of the photons. The cross section thus obtained would receive most of the 
contribution from the region in phase space where the Ig^^l are small (this also corresponds to 
small scattreing angles for the proton) because of the 1 / factors from the photon propagator. 
Note that, there is a kinematic lower bound on |g| 



q^ < 



where E is the energy of the proton in the center of mass frame and rrip its mass (see Appendix 
C for the derivation). 

The fact that most of the contribution to the cross section comes from the small \qf\ region 
means that we can evaluate the contribution to the amplitude from the 77 — )■ X part of the 
diagram in the \qf\ — )■ limit. This is the so-called Equivalent Photon Approximation (EPA). 
This amounts to treating the photons as real with only transverse polarizations while doing 
the the 77 —t- X part of the calculation so that the total cross section can be written in the 
factorized form, 

(35) 



a 



r-2 



2Mz 



Here W = a/I is the photon-photon center of mass energy and a^^ is the photon level cross 
section. Sq^j^, the survival probability for diffractive photon exchange processes, is the prob- 
ability that the proton remains intact and is not broken due to subsequent inelastic QCD 

0.9 following the theoretical calculation in Ref. 



interactions. We take Sq^^ 



11 



The func- 



tion dL/dW contains all the details of the proton electromagnetic form factors and also the 
integral over 1/qf factors of the photon propagators. A detailed calculation of dL/dW using 
EPA leads to the following expressions (see Appendix D in Ref. [19]), 



dL 
dW 

C 







aE 
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2Wfix)f 
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q'Gi,)l{Aml - q' 



(36) 
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Here x = ^ and s = AE"^. While the lower limit of the integration is set by kinematics (see 
34) we take the upper limit to be g^^^. = 2 GeV. Beyond = 2 GeV, the form factors 



eq 



in eq. 33 become very small so that the contribution to the integral is negligible. 



To understand the physical meaning of dL/dW we can multiply both sides of eq. 35 by Cp, 
the proton luminosity. Then we find that the luminosity function is the ratio of the differential 
photon luminosity dC-y/dW and the proton luminosity, 

dL/dW = 5*^. (37) 

Note that here L is unitless and C-y^p has the usual units m~^s~^. We plot the photon lumi- 
nosity function in Fig. [3j We find that, 

„ r'" dL 



For a particular process 77 — )■ X this number gives an upper bound on the ratio, 

(j(77 — )■ X) 

if cr^-y is a constant or decreasing function of W as is the case usually for SM processes. Thus 



from a knowledge of (T-y-y(77 — )■ X) one can estimate a (pp(77 — )■ X)pp) using eq. 38 

6 Theoretical cross sections 

In this section we present the cross section for pp(77 — )■ 77)pp and pp(77 — )■ ZZ)pp (see 
Fig. [T| with the proton-proton center of mass energy equal to 14 TeV. We will consider only 



the light higgs case in Section 2.1 taking all fej/A^ = 1 TeV ^ and all q/A^ = 1 TeV ^ in 
eqs. |9]and|2} With these values for the couplings and using eq. 30 and [31] for / = 0, we get 
the unitarity bound Kffj^ = 1220 GeV and = 1260 GeV respectively. 



We have evaluated the cross section with and without the form factor in eq. 28 For the 



form factor we have taken Aj^^^^ = , m = 2 and n = 1. We have taken the Higgs 

mass rrih = 120 GeV. We have evaluated the cross section with the cut W < W^ax and varied 
Wmax- This cut is important as the ambiguities due to the choice of form factor become more 
important for large values of Wmax- For the ppi'yy — ?■ '^l)pp process we show the results 
with and without the form factor in Fig. |4j We have checked that the contribution from the 
77 —;■/;.* —> 77 process due to the presence of the operator Oq is small compared to the total 
cross-section as is expected from the arguments in Section |4j 
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Figure 4: The ppi'^'j — )■ 'y'j)pp cross section we obtain as a function of Wmax with and without 
a form factor. For the form factor we use in eq. 28 with m = 2, n = 1 and Ay^^^ = K^J^^^ . 
We have taken hi/K^ = 1 TeV'^ and c^/A^ = 1 TeV-Mn eq. Isl the Higgs mass rrih = 120 GeV 



and the proton-proton center of mass energy equal to 14 TeV. 

We show the results for the pp(77 — t- ZZ)pp process with and without the form factor 
in Fig. [sj In Fig. |6]^left) we show the 77 — > ZZ cross section we obtain as a function of 
Wmax without any form factor. We show separately in the same figure the contribution 
due to the operator Oq through the process 'j'j ^ h* ZZ. We see that the non-local 
contribution due to Oq dominates at low energies. The contribution to the cross section due 
to the other operators, however, grows more rapidly with s (as a^^ ~ s^) compared to the 
Oq contribution. At higher energies Of ^ contributes most to the cross-section. As mentioned 
earlier the coupling oq can be accurately measured by measuring the /;.—)■ 77 partial width so 
any deviation would indicate the presence of higher dimensional operators. We also show in 
Fig. [6]^left), the 77 — )■ Z^Z^ and 77 — )■ Z^Zt contributions to the cross section. As explained 
before, for longitudinal Z production, the main contribution comes from the operators Oq, C^f f 
with Oq contributing dominantly at low energies and Of ^ contributing dominantly at higher 
energies. For transverse Z production only the operators Oq, C^^s g contribute significantly 
with the dominant contribution coming from Of ^. 
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Figure 5: The ppi'j'j — > ZZ)pp cross section we obtain as a function of W^ax with and without 
a form factor. For the the form factor we use in eq. 28 with m = 2, n = 1 and = A^^^^. 

We have taken bi/A"^ = 1 TeV'^ and c^/A^ = 1 TeV^ in eq. Isl the Higgs mass rrih = 120 GeV 



and the proton-proton center of mass energy equal to 14 TeV. 

In Fig. |6](right) we have the signal cross section curves as in Fig. [6|^right) but with form 
factors. The contribution due to Oq is shown without any form factor suppression. This is 
because we want to show the pure contribution of the operator Oq so that any deviation can 
be interpreted as the presence of higher order effects in s/A^ (as explained in Section |4| using 
a form factor would amount to assuming such higher order corrections). 

Finally let us state how the contributions of the different jjZZ operators in eq. [7] can, in 
principle, be resolved. We have already seen how looking at the final polarization of the Zs 
can be used to distinguish the contribution of Of ^ and Of ^ from the other 'j'jZZ operators. 
Another fact that can be used is that only for the operators Of ^ and Of ^ are the amplitudes 
spherically symmetric. Thus the operators in eq. [7] can be divided into four categories: 

those that contribute mainly to the Z^Zl mode and give spherically symmetric amplitudes 
(only Of^), those that contribute mainly to the Z^Zl mode but do not give spherically 
symmetric amplitudes (only Of^), those that contribute mainly to the ZtZt mode and give 
spherically symmetric amplitudes (only Of^), and those that contribute mainly to the ZtZt 
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Figure 6: The pp{'yy — )■ ZZ)pp cross section we obtain as a function of Wmax without any 
form factor. We have taken bi/A'^ = 1 TeV~^ and q/A'^ = 1 TeV~'*in eq. [sj the Higgs mass 
rrifi = 120 GeV and the proton-proton center of mass energy equal to 14 TeV. We show 
the total cross section (solid), the ZtZt production cross section (big dashes), the Z^Z^ 
production cross section (small dashes) and the total cross section only due to the operator 
Oq (dotted) through the 77 —)■/;.* —j- ZZ process. 

mode but do not give spherically symmetric amplitudes ( O^^ — O^^). Note that resolving the 
contributions of the different operators would require higher luminosity than just detecting 
the presence of QNGCs, but we will not go into the experimental feasibility of such studies. 

7 LHC signal search strategy 

As explained in Section [5] the final state protons in diffractive processes are scattered at small 
angles. To detect such protons very forward detectors have been proposed both for the ATLAS 
and CMS detectors (see Ref j4|j). It has been proposed that such detectors should be placed 
at distances of 220 m and 420 m from the interaction point where the distance is along the 
circular beam line. To give an idea about these distances, at 220 m the beam line curves away 
from the tangential direction at the interaction point by about 6 meters. The protons lose 
a small fraction of their energy in the diffractive process and experience a small deflection 
from the beam axis. As the deflection is very small, the LHC magnets continue to curve the 
protons along the beam pipe but they do move away form the beam axis and out of the beam 
envelope because of the deflection. Thus detectors close to the beam axis (a few millimeters 
away) would be able to detect the protons. It is also important to note that particles other 
than protons would never be detected in these detectors as they have a different cyclotron 
radius. Thus these detectors effectively use the LHC magnets as a spectrometer. 
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As the detectors need to be close to the beam axis radiation hardness is a requirement that 
the detector must fulfill. This along with resolution requirements makes 3D silicon detectors 
ideal as proton detectors. From the measurement of the position and track direction at the 
detectors the momentum four vector of the proton can be reconstructed by inverting the 
transport of the proton due to the LHC magnet optics. Thus it is possible to measure the 
fraction of energy lost by each proton, C,i and thus measure the invariant mass of the central 
system (also called the proton missing mass), 

W = (39) 

where s = (14 TeV)^ for the 14 TeV LHC. The 220 m detectors detect protons with smaller 
deflection, and thus smaller ^, than the 220 m detectors. As higher invariant masses would 
correspond to higher ^s, the 420 m detector is sensitive in the low mass region whereas the 
220 m detector is sensitive in the high mass region. The 220 m detector is thus crucial for the 
kind of study we are doing in this paper where most of the signal contribution comes from 
events with high W. Mass resolution between 2 GeV and 3 GeV for low energies and about 



5-6 GeV for the highest photon energies can be achieved by these detectors 20 . 

As only events with two intact protons are accepted, the only background processes can 
be those in which the proton emits a light particle with no electromagnetic or color charge 
and remains intact. Thus the proton can emit a photon or a pomeron. In pomeron fusion 
processes, also called 'Double Pomeron Exchange' (DPE) processes, the pomeron, in general, 
breaks into fragments. Thus whenever we would write down a process pp[DPE — )■ X)pp, X 
being a particular final state, it would be implicit that this is an inclusive process where other 
particles (pomeron fragments) are also present. These detectors can be used to test if an event 
is exclusive or inclusive, where by an exclusive event we mean an event where no other particle 
in addition to the final state particles is produced. This can be done by matching the invariant 



mass measured by the proton detectors (using eq. 39 ) with the invariant mass measured by the 
central detectors. Also, in exclusive events the px of all the final state particles excluding the 
protons (which carry very little Pt) must add up to nearly zero. Thus if only such exclusive 
events are accepted the underlying process can only be a an exclusive pomeron fusion process, 
usually called 'Central Exclusive Production' (CEP) and the inclusive DPE background can 



be reduced. Including both the 220 m and 420 m the acceptance range for is [20] , 

0.0015 < < 0.15 (40) 



Using eq. |40] we see that only events with 21 GeV< W < 2100 GeV are accepted by the 
detector. 
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A potentially important background contribution is from overlap events. If the signal 
event is pp — )■ pXp, an overlap event would be defined as the coincidence of an event where 
the central system X is produced with one or more diffractive events in the same bunch 
crossing. Processes like [p] [Xp], where one of the protons is produced in an interaction different 
from the central process (the square brackets grouping the final particles produced in the 
same scattering process) or processes like and [p][X][p] where both the protons are 

produced in an interaction different from the central system, can fake signal events. To reject 
such background events the forward detectors would be equipped with timing detectors which 
would have a resolution of the order of 10 ps |4| (note that the LHC bunch length is about 
250 ps at 1-sigma). These detectors would be able to reconstruct the vertex position of the 
two protons assuming they are produced at the same interaction point. This vertex would 
not match the vertex for the central system, X, for the fake events and thus such background 
events can be rejected. Another way to reduce the overlap background is by matching the 
net invariant mass of the central system measured by the central detectors with the values 
obtained by the forward detectors. We have already discussed how this can be done in the 
context of testing whether an event is exclusive or not. Overlap backgrounds are of great 
importance when the inelastic production cross-section for the central system X is large as 
is the case for the dijet background to diffractive H ^ bb production [2l] (where the jets are 
misidentified as 6-jets) but is of much lesser importance in our case. Let us now discuss the 
search strategy for the pp(77 — ZZ)pp and the ppipf^ — )■ ll)pp processes. 

7.1 ppijj ZZ)pp process 

As we do not perform a detailed detector simulation and wish to make only an estimate of the 
detector level cross-section, we will look at the p + {ZZ — )■ 4/) + p final state (/ = e, /i) that 
is most free from experimental complications. Other final states involving hadronic Z decay 
modes may well turn out to be more sensitive to our observable but ascertaining this would 
require a more rigorous experimental analysis. 

The most important background process is pp{DPE — )■ ZZ)pp in the SM. As mentioned 
before this is an inclusive process where the final state has pomeron remnants in addition 
to the Z pair. The DPE background has been computed using the Forward Physics Monte 
Carlo (FPMC) |23| which uses the Ingelman-Schlein (IS) model j22] for inclusive diffraction. 
In this model the cross-section of a process like pp{DPE — )■ X)pp is computed by convoluting 
the cross-section of the partonic hard process ij — )■ X being the partons) with diffractive 
parton density functions (DPDF) measured at HERA. The DPDF describes the probability of 
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extracting a certain parton with a given longitudinal momentum fraction from the proton. The 
DPDF itself can be expressed as a product of the pomeron flux, measured in other diffractive 
processes, and the probability of extraction of the parton from the pomeron which breaks 
into fragments. Whereas the IS model describes diffractive data at the Hadron-Electron Ring 
Accelerator (HERA) very well its theory prediction for diffractive dijet production at the 



Tevatron is larger by a factor of 10 24 . This is usually attributed to the fact that there is 
some probability for the protons in a DPE process to have a subsequent inelastic interaction 
which breaks the proton. Thus the cross-section computed in this model must finally be 
multiplied by the survival probability, Sj^p^, which is the probability that there is no further 
inelastic interaction between the protons. The maximum value of the survival probability 



reported in the literature is about Sj^p^ = 0.06 25 . 

For the pp{DPE ZZ)pp process that we are interested in, there are two possible 
partonic sub-processes, the qq — ?■ ZZ sub-process via t-channel quark exchange of a quark 
q and the gg — ?■ ZZ sub-process that is induced by fermion loop^ The quark component 
of the DPDF gives the dominant contribution in processes involving diffractive vector boson 



production like the one we are interested in 26 . We apply the following cuts to the signal 
and background |8 ,20 , 



0.0015 < < 0.15 (41) 
W > 300 GeV. (42) 

The first cut above is just the ^-acceptance cut for the forward detectors, and the second cut 
has been applied mainly to suppress the SM loop background discussed later. As explained 
in Section |6| most of the contribution to the signal cross-section comes from high energies so 
that the second cut hardly affects the signal. With these cuts the cross-section we thus obtain 
for pp{DPE — )■ ZZ)pp from FPMC including all Z decay modes is 1.4 fb. Applying the above 
cut the signal cross-section (with form factor) is reduced from 8 fb to 3 fb in the light Higgs 
case for fej/A^ = 1 TeV~^ and Cj/A^ = 1 TeV~^, where the ^-acceptance cut is responsible for 
most of the reduction. 

As the DPE background discussed above is inclusive unlike the signal it can be further 
reduced by testing if the events are exclusive. This can be done by matching the four-lepton 
invariant mass measured by the central detector with the invariant mass measured using eq[39j 
Also the net px of the four leptons must add up to nearly zero (as the protons carry hardly 

^For the computaion of the partonic cross-section of the gg — ?> ZZ sub-process, the gg ^ h* ^ ZZ 
process has not been included. This contribution is known to interfere negatively and decrease the total cross- 



section 27 so the background cross-section would have been lower had this contribution been incorporated. 
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any transverse momentum) for exclusive events. 

Now let us consider the background contribution from the SM loop process pp(77 — )■ 
ZZ)pp. The cross section (7^^(77 — )■ ZZ) has been evaluated in Ref. 28 to be roughly- 



constant, around 300 fb in the range 300 GeV < r < 2100 GeV, where the upper limit is 



equal to the upper limit obtained by applying the ^-acceptance cut in eq. |46j Thus we find 

for the proton level cross section, 

/■sioo /.2100 
a = SIed I -^cr^^dT ^ 0.9 x 300 / —dr ^ 0.1 fb, (43) 
J2MZ dr J2MZ dr 

where we have imposed the cut 0.0015 < < 0.15. We find therefore that this background is 
negligible compared to the signal. 

So far in all the instances where we have considered signal or background cross-sections we 
have not taken into account the branching ratio of the Z to leptons and detector efficiencies. 
Taking the lepton identification rate to be 90% |j29j and the proton detection efficiency in 
the forward detectors to be 85% |4| we obtain for the effective cross-section we expect the 
detectors to measure, 

a,ff = 0.56 BiZ ^ llfath, (44) 

where ath is the theoretical cross-section including all Z decay modes. 

Finally, let us mention a possible complication that may arise because of the fact that the 
high energy Zs we are considering would be boosted in the lab frame. This would cause the 
leptons to be coUimated along the direction of motion of the parent Z. This may give rise 
to complications in detection of some electron pairs for which the two electrons are not well 
separated from each other (there is no such issue with muonic decays as muon separation is 
always efficient for the energies we are considering). We will not try to estimate this effect 



(see for example for a 30 more detailed discussion) but in our estimates of sensitivity in the 
next section, to give a conservative estimate, we will provide results considering only muonic 
decay of the Zs in addition to results considering decay of the Zs into both electrons and 
muons. 



7.2 ppi'jj — > ll)pp process 

For the ppi'j'j — )■ 'y'y)pp process we require the presence of two photons and two protons in 
the final state. Again the main background is from the pp{DPE — ?■ '~f'y)pp process. The 
pp{DPE — '~f'~f)pp cross-section can be estimated from the pp{DPE — )■ ZZ)pp cross-section 
by using the fact that the t-channel quark exchange sub-process qq ZZ/'-fj is the dominant 
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partonic sub-process 26 . As the diagrams for this partonic sub-processes in both the cases 



are the same except for the two outgoing vertices and external legs in the limit of massless 
quarks we get, 



a{pp{DPE — i- '~f'~f)pp) 
a{pp{DPE ZZ)pp) 



[l-AMl/s)- 



-l/2_ 



(fl'/2c^ 



(45) 



where Vq and are the vector-like and axial vector-like couplings of the quarks to the Z boson 
and Q is their electric charge. The sum is over the three light quarks for which the probability 
of diffractive extraction from the proton is significant and we assume that the diffractive PDFs 
for three light quarks are equal. The kinematic factor on the RHS, which is almost unity at 
high energies, arises because the Z-boson unlike the photon is massive. Taking Vd^s = —0.35, 
Vu = 0.20, au = 1/2, ad = —1/2 and s = (500 GeV)^, we find this ratio to be 0.3. We apply 
the same cuts as in the case of the pp{'~f'~f ZZ)pp process, that is. 



0.0015 <^i< 0.15 
W > 300 GeV 



(46) 
(47) 



and using the above ratio we obtain this background cross-section to be 0.4 fb. Again this 
inclusive DPE background can be further reduced by requiring the two photon invariant mass 
to match the missing mass evaluated using eq. 39 and by demanding that px-yi = —PT-y2 within 
experimental resolution. 

The SM loop induced ppi'yy — )■ 77) process in this case has a cross-section that is C(O.Ol) 
fb and can be ignored [3l||32] . An experimental background contribution can come from mis- 
identification of jets as photons in the pp(77 — )■ jj)pp process. The total inclusive diffractive 
dijet cross-section at the LHC has been computed by the diffractive Monte Carlo generator 



DPEMC 33 to be 4x10^ fb 34 in the IS model with the cut Et > 25 GeV for the jets. 



Taking the rejection factor of jets for photon identification to be 5000 29 we get a background 



cross section of about 2 fb which is already smaller than the signal cross-section (with form 
factor) of 17 fb in the light higgs case for bi/A"^ = 1 TeV~^ and q/A^ = 1 TeV~^. Further 



cuts like the ly-cut in eq. 47 and requiring px-yi = — PT72 within experimental accuracy should 
completely remove this background. 

The effective detector level cross-section is again smaller than the values mentioned so far. 
Taking the photon identification rate to be 90% [29) and proton detection efficiency in the 
forward detector to be 85% |4] we get. 
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0.69 ath 



(4J 



25 



Couplings 


Process 


InteEfrated 
Tjiiminositvf fh~^l 


^ ' oos 






Confidence 
Tjevel I siffma 1 


^^ase i. [oov yjiGv ) 

^dbfi i. ^^i.o lev ) 


77 — ^ 77 


1 


IZ. 1 




n Q 
u.o 


\ 1 n 
> lU 


77 — ^ 77 


'inn 
ouu 


1 QQ 1 
loo. 1 




89 8 


o.z 


Uase i. [odO Lrev ) 


77 — y Z Z 


Qnn 
oUU 


v /I /I 
( .4(^i.y ) 


1 
i 


1 ^'n Q^ 
1(0.3) 


A Q ('9 1 ^ 

4.o(^/. i ) 


Uase i. [lou Lre V j 
uase i.(^ouu vjrev j 


77 Z Z 


oUU 


11/1 /'o Q^ 
ii.4(^/.o ) 


i 


1(0.3) 




77 — ^ .Z'.Z' 


Qnn 
oUU 


4D. 0(^11. ( ) 


1 
1 


1 (V\ Q\ 

1(0.3) 


\i n('8 1 ^ 
>lU(^o.l J 


uase z.(^(UU vorev ) 


77 — ^ .Z'.Z' 


Qnn 


1/1 S/'Q 7^ 

14.0(^0. ( ) 


9 


1 ^'n K\ 
1(0.5) 


o.oyo. 1 J 




77 — ^ ^ Z 


Qnn 
oUu 


c;i qH 8^1 
01. 0(^1/. ) 




z 


1 ^'n 




Pn«P — 1 n TpV 
l^aoc O. IVX — -L-'J lev 


y y ' y y 


1 
1 


1 'i ^ 
10.0 




n '\ 

u.o 


^1 n 

^ lU 


Case 3: At = 2.4 TeV 


77 — )■ 77 


300 


118.2 




82.8 


3.9 


Case 3: At = 900 GeV 


77 — )■ ZZ 


300 


12.6(3.2) 


1 


1(0.3) 


6.4(3.6) 


Case 3: At = 700 GeV 


77 — )■ 


300 


39.6(9.9) 


1 


1(0.3) 


>10(7.1) 


Case 4: (1.9 TeV)-^ 


77 — > 


300 


5.3(1.3) 


1 


1(0.3) 


3.3(2.1) 


Case 4: (2.2 TeV)-^ 


77 — 7- ZZ 


300 


3.9(1.0) 


1 


1(0.3) 


2.2(1.1) 



Table 1: The expected number of observed events Nobs and the signal significance for both 
the processes for different integrated luminosities. The expected number of observed events 
is evaluated using Nobs = crlJJ^"'''L where L is the integrated luminosity, and to evaluate the 
signal significance the background is assumed to follow a Poisson distribution with mean 



Nb — cr'y/fL. The signal contribution has been evaluated with a form factor as in eq. |28| taking 



m = 2, n = 1 and Aj"'^^^ = A^^^^ . For the 77 — )■ ZZ process the values in the parentheses 
show the results if only muonic decays of Z are considered. 



Note that for both the pp(77 — )■ 77)pp and pp(77 — )■ 77) ZZ processes we have ignored 
above the effects of the basic detector acceptance cuts pr > 10 GeV and r] < 2.5 for the 
leptons and photons. As the dominant contribution to the signal cross-section is central and 
from high energies, these cuts are expected to have a very small effect. 



8 LHC sensitivity to QNGCs 

Using the LHC search strategy for ^^(77 — )■ ■yj)pp and ppipf^ — )■ 'y'y)ZZ signals outlined in 
the previous section we can now report the expected sensitivity of diffractive photon fusion 
at LHC to QNGCs. Table [T] shows the expected number of observed events Nobs and the 
signal significance for both the processes with different integrated luminosities. The expected 
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number of signal, background and observed events are evaluated using, 

Ns = aZ'r^Ci 



Nobs = Ns + Nb (49) 



where Cint is the integrated luminosity and aeff is the effective cross-section defined by eq. 44 



and eq. 48, after taking into account detector efficiencies. The signal contribution has been 
evaluated with a form factor as in eq. 28 taking m = 2, n = 1 and Aj^^^^ = A]^^^^ . We 
can simply add the signal and background events to get the total number of events expected 
to be observed in eq. 49 because the interference with the background is very small. The 
interference with the DPE background is small because the interference between DPE and 
photon exchange diffractive processes is in general small and the interference with the SM loop 
background is small because unlike the signal this background gets most of the contribution 
from the \ow-W region. In order to quantify the signal significance we evaluate the probability, 
a, that the background has not fluctuated to give a number of events greater than or equal 
to Nobs assuming that it follows a Poisson distribution with A*";, as its mean. The confidence 
level expressed as a particular number of sigma deviations is given by \/2 erf~^(a;) where erf() 
is the error function. We find the sensitivity for four different physically interesting ways of 
choosing the relative value of the QNGCs. 



CASE I: bi/A^ = l/A\ a/A^ = l/A^ 



We find that the 77 production process is by far the more promising of the two processes 
for probing QNGCs. As we can see from Table [T] even with integrated luminosities as low as 
Cint = 1 fb"^, couplings as small as 1/(850 GeV)^ can be probed with large significance. With 
high integrated luminosity (300 fb~^) couplings as small as 1/(1.8 TeV)^ can be detected with 
more than 5 sigma significance. There are possible cuts that can remove inclusive events as 
discussed in the previous section, which may substantially reduce the DPE background. If 
this is possible the 77 — )■ 77 process can be sensitive to even smaller couplings. Note that 
a coupling with value 1/(1.8 TeV)^ does not necessarily mean that the energy scale of new 
physics is 1.8 TeV. If dimensionless couplings less than unity or loop factors are present, for 
instance, the scale of new physics would be lower. 

The ZZ production process requires very high integrated luminosity. For this process we 
give in addition to the results assuming Z decays to both electrons and muons, the results 
considering only the muonic decays in parentheses. For Cint = 300 fb~^ as one can see from 
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Table [T] the smallest couplings that can be detected with more than 95 % confidence level are 
about 1/(850 GeV)^. If we require the detection of at least 10 signal events these values are 
1/(750 GeV)^ and 1/(500 GeV)^ considering respectively decays to both electrons and muons 
and only muonic decays. 

CASE 2: Resolving QNGC contributions to ZZ production from the contribution due to Oq 



As mentioned before the dimension-6 operator Oq contributes to the signal through the 
77 — )■ /i* — )■ ZZ process. The value of this coupling can be obtained from the /i — )■ 77 
partial width measurement. For this case we consider this contribution to be part of the back- 
ground and take all the 6i/A^ = 1/(850 GeV)'*and all the Cj equal. We then try to find the 
smallest QNGC couplings q that can be detected. As we want to separate the Oq contribution 
from higher dimensional contributions, we do not use any form factor for the evaluation of 
the 77 —;■/;.* ZZ cross-section due to this operator as using the form factor is equivalent to 
including higher dimensional corrections (see Section |4]). We find that the smallest couplings 
that can be detected for this case with 300 fb'^ data to be 1/(700 GeV)^ (1/(500 GeV)^) 
considering Z-decays to both electrons and muons (to only muons). 



CASE 3: Graviton exchange in extra-dimensional model 



For this case we assume that the QNGCs arise from the effective dimension-8 operator due to 



virtual graviton exchange in extra-dimensional theories described in Section 2.3 The relative 



couplings of the QNGCs are thus fixed by expanding the operator in eq. 16 and the only 
adjustable parameter is Ay. 

As one can see from Table [T] we find that the 77 — )■ 77 process can detect this operator 
with only 1 fb~^ data for A^ = 1.0 TeV. For high luminosities (300 fb~^) the maximum value 
of At that can be probed by this process in the DPE background is about At = 2.4 TeV. 



Note that our results differ from and are less optimistic than the results of Sahin et al 32 



who do not consider the DPE background and more importantly use a far less restrictive ^- 



acceptance cut. As explained after eq. 40 because of the ^ < 0.15 acceptance cut only events 



with W < 2.1 TeV are accepted. In Ref. 32 , on the other hand, events with ^ as high as 0.5 
are accepted which corresponds to W as high as 7 TeV and most of the contribution to their 
signal comes from the high 1^ events; protons with C, > 0.15 can, however, not be detected by 
the forward detectors [8,20 . For the 77 — )■ ZZ process At as high as 900 GeV (700 GeV) 
can be probed with 300 fb~^ data considering Z-decays to both electrons and muons (to only 
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muons). As already mentioned in Section 2.3 the gg'j'j/ggZZ operators {g being a gluon) 



that arise from expanding the operator in eq. 16 are expected to give a contribution to the 



exclusive ppjj/ppZZ final states via central exclusive pomeron fusion but this contribution 
is expected to be negligible relative to the diffractive photon fusion contribution. 

Establishing the presence of QNGCs would give very important complementary evidence 
for virtual graviton exchange because it is possible in this case to uniquely trace back to the 
underlying dimension 8 operator involved. Our final sensitivity results show, however, that for 
the particular diffractive processes we have studied for probing QNGCs the largest At that 
can possibly be probed (2.4 TeV) has already been ruled out by dijet constraints from the 
36 pb~^ CMS data in Ref. |l4j where the constraint Aj- > 3.8 TeV has been derived. Thus 
diffractive photon fusion will not be able to probe A7- values still allowed by experimental data. 

CASE 4: Higgsless case 



As we discussed in Section 2^ in the higgsless case we expect only the following two op- 
erators to be important, 

^qTgc'' = ^-^^^^^F.^F'^'Z^Z^ + ^^^^^^^^F,^F^^Z,Zr (50) 



We take a^' = in eq. 15 and find that the pp{'~f'~f — )■ ZZ)pp process is sensitive up to 
couplings as small as 1/(1.9 TeV)^ if we require more than 95% confidence level, a huge 
improvement over existing limits. Our sensitivity estimates agree well with those obtained by 
Royon et al. in Ref. [s], once we translate to their convention for parametrization of these 
couplings. 

Higgsless models are usually associated with strong electroweak symmetry breaking (EWSB) 
scenarios. The operator coefficients in such theories can be estimated by Naive Dimensional 



Analysis (NDA) (see Refs |35||36|). In our case, using NDA, we find /A^ = / {IQi:'^ Kf) , A, 
being the scale of the strongly coupled sector]^ The above mentioned estimates tell us that 
our process is sensitive to A^ < 100 GeV. As a strong sector at such low energies is already 
ruled out by experiments, our process, unfortunately, cannot probe realistic scales for strong 
EWSB. 



^Note that the ZZZZ couphng which appears at the dimension 4 level in the chiral lagrangian, from 
operators like c {Ty[{D ^X)^ D^'Y.])'^ , is less suppressed (c ~ v"^ / « l/(167r2)) than other QNGCs. 
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9 Conclusions 



We have listed all possible operators contributing, at the lowest order, to Quartic Neutral 
Gauge Couplings, quartic gauge couplings involving only the photon and the Z boson and have 
studied the sensitivity of measurement of these couplings in diffractive photon fusion processes 
at the LHC. These couplings are interesting because the lowest order contribution they receive 
is from dimension 8 operators in scenarios with a light Higgs and, in higgsless scenarios, from 
dimension 6 operators (with the exception of the ZZZZ coupling which receives dimension 4 
contributions in this scenario but we have not focussed on this coupling in this work in any 
case). Thus new physics processes which do not contribute through operators of the lowest 
possible dimension can be probed by measuring these couplings. One specific example that 
we have considered is virtual graviton exchange in extra dimensional theories where the lowest 
dimension operators generated are of dimension 8, and these include operators contributing 
to QNGCs. 

Thus measurement of QNGCs in any experimental process would be interesting, but 
in this work we have studied their measurement in diffractive photon fusion processes like 
PPill ~^ ll)pp si-nd pp{'y'y — )■ ZZ)pp. The protons in these processes remain intact and 
scatter diffractively with very small scattering angles. These can be detected by very forward 
proton detectors that have been proposed for both the ATLAS and CMS experiments. As we 
argue the detection of the two 'j/Zs in the central detectors along with the detection of the 
protons in these forward detectors would indicate the existence of QNGCs like the 7777 and 
^'yZZ couplings, as this is the only feasible new physics possibility that can lead to such a 
final state. The only other possibility is pp{CEP ZZ)pp, where CEP stands for Central 
Exclusive Production, is a process that takes place when pomerons fuse exclusively (that is 
without breaking into fragments) to give the ZZ/77 final state. Such processes are, however, 
expected to have a much smaller cross-section when compared to photon fusion processes. To 
calculate the cross-section for the pp(77 — )■ 77)pp and pp(77 — >■ ZZ)pp processes we convolute 
the cross-section of the 77 — )■ 'j^/ZZ sub-process with the 77 luminosity function obtained 
using the Equivalent Photon Approximation. The amplitude of the 77 — )■ 77/ZZ sub-process 
grows with energy because of the non-renormalizable couplings involved and we unitarize this 
using appropriate form factors. We have argued that our final sensitivity results for A will 
not change much for a different choice of form factor than ours. 

Before we summarize our results on the sensitivities, note that QNGCs are very weakly 
constrained by existing data. Whereas no constraints exist on 7777 couplings, the 77ZZ 
couphngs are constrained by direct search results from LEP to be smaller than about 1/(100 
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GeV)^ (1/(50 GeV)^) and by precision measurements to be smaller than about 1/(270 GeV)*^ 
(1/(420 GeV)^) in the light Higgs (higgsless) case. We have found in this study that diffractive 
photon fusion at LHC can improve these sensitivities by many orders of magnitude for the 
coupling, and can probe couplings as small as 1/(850 GeV)^ (1/(1.9 TeV)^) with 300 
fb^^ integrated luminosity for the light Higgs case (higgsless case). We find, however, using 
an NDA estimate, that the values in the higgsless case correspond to a scale lower than 100 
GeV for the strong sector which is already excluded by experiments. The 7777 coupling can 
be probed even more sensitively and values as small as 1/ (1.8 TeV)^ can be measured with the 
same integrated luminosity for the light Higgs case. For the specific case of virtual graviton 
exchange in theories with large extra dimensions we find that the highest scale that can be 
possibly probed (about = 2.4 TeV by the pp(77 — > 77)pp process with 300 fb~^ data) 
has, unfortunately, already been ruled out by the latest constraint from CMS dijet data which 
puts the bound At > 3.8 TeV. 
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Appendix A: 777^, jZZZ and ZZZZ couplings 

We first consider the operators in the light higgs case in eq. |2] when written in terms of 
the fields A and Z give rise to 7777, j'jZZ, j'j'-fZ, 'jZZZ and ZZZZ couplings. We already 
wrote the Lagrangian for the 7777 and 'j'jZZ couplings in eq. [7} In this Appendix we will 
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present the 'j'j'jZ, •yZZZ and ZZZZ couplings. Using eq. |6]we get from eq. |2} 

-y-y-yZ 777^ 
^yyjZ _ "1 p p^iu p ypu I ^ p pup p yuv 

-tZZZ ^2 ^"/ZZZ ^-/ZZZ 

Clfa'c = ''-^^fF.^Z^^Z.Z^ +''-^^^F,^Z^^Z,Z^ +''-^F 

■yZZZ 

^'^^F^.Z^^PZ^^Z'^- 

ZZZZ T\,fi ZZZZ ZZZZ j\yT2 

Z Z Z nZ z z z 

n ^^4 ^i^y^ ^pa^ ~\ ^^4 ^P^y^ ^pa^ y->^) 

where for the 777Z couphngs we get, 

aV^^ = -4s^C^C8 + 4S^C^(C9 + Cio) + (2s^„C^ - 2sj^Ct„)(cii + C12) 

al'^'^^ = -4s^c^ci3 + 4s^c^(ci4 + C15) + (2s^c^ - 2s^Cu,)(ci6 + C17), (52) 
for the 'jZZZ couphngs we get, 

aj = 2SwCwC3 + 2s^C^C4 (c^ ^w)^^ 

^ZZZ r) I O 

al^^^ = — 4s^c^t,C8 + 4s^c^(c9 + Clo) + (2s^Cu, — 2s^c^)(cii + C12) 

aj^^^ = -4s^c^ci3 + 4s^c^(ci4 + C15) + (2sj^c^ - 2s^4)(ci6 + C17). (53) 

and for the ZZZZ couphngs we get, 

af^^^ = C1 + C2 

qZZZZ _ ^^Q^ -)- + C^SttjCs 

ZZZZ _ „2 , 2 

flg — s^ce + c^cr 

af^^^ = stcs + ct(c9 + cio) + c^s^(cii + C12) 

a™ = stci3 + ct(ci4 + ci5) + 44(ci6 + ci7). (54) 



As explained in Section 2.2, any U(l)em invariant operator, constructed using and F^y 
fields, is an allowed operator in the higgsless case. Thus, for the higgsless case we will get 
the same operators as above but now Cf would arise from dimension 4 operators while 
^7zzz^ Q-iZzz ^ qzzzz (jzzzz -^Q^ifj arise from dimension 6 operators. 
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Appendix B: Derivation of unitarity relation 



In this appendix we derive the expression for the unitarity bound for the processes 77 — )■ 77 
and 77 — )■ ZZ. Applying the optical theorem to the 77 — )• 77 process tells us, 

Im(Al (7172 ^ 7172)) , . 
= (7(7172 — )■ everything) 

s 

= cr(7i72 7(61)7(62)) + ^ ct(7i72 Z{e3)Z{e4)) 

+A (55) 



where 7, denotes 7(^1, €«) and A is a positive number that accounts for all the other con- 
tributions to the RHS of eq. 55 and the cross section for the 77 — )■ VV process is given 
by, 

/3v 



>w 



J dncu \M{^il2-^VV)\\ 



The amplitude can be expanded into partial waves as follows, 



(56) 



(7172 ^ 7172) = 167r^(2J + l)6jPj(cos0) 
A^(7i72^^^) = 167r^(2J + l)ajPj(cos0). 



where /3- 



V 



1 



4Mt 



. For the forward scattering in the LHS of eq. 
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(57) 



we must put = 0. 



Using eq. 55, 57 and 56 and the following property of Legendre polynomials, 

2 



I 2n + 1 



gives. 



( Im(60)' - Im(6z) + J^iMbi))^ + PwYl l^'l' + = 0- 



(58) 



(59) 



The first two terms in eq. 59 should be evaluated taking the initial polarizations to be exactly 
same as the final polarizations, and Si is the positive contribution from every other source. 



Eq. 59 is a quadratic equation for lm{bi). The equation must have real roots and thus must 



have a positive discriminant. This gives the condition, 

(Re(60)' + /3$^|a/P + 5/< ^. 



(60) 
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Appendix C: Kinematic hound on photon virtuality 



First let us derive the kinematical limits on qf. From conservation laws we must have 
q = {E — El,p — ff^). Substituting = ^jE^l^^rrfi ^ m being the mass of the particle emitting 
the photon, we obtain for m <^ E^, 

q^ = -4EE[sin'^^-'^cos9,. (61) 

Here 6i is the angle between p and p^. In the expression above the first term dominates. 
As most of the contribution to the amplitude comes from the small Ig^^l region, ignoring the 
second term above we see that we must have small 6. We thus obtain the following kinematical 
bound on qf, 

* < -m^y '''' 
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